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Abstract 



We place bounds on long-lived primordial relics using measurements of the diffuse gamma 
ray spectrum from EGRET and COMPTEL. Bounds are derived for both radiative and 
hadronic decays with stronger bounds applying for the latter decay mode. We present an 
exclusion plot in the relic density-lifetime plane that shows nontrivial dependence on the 
mass of the relic. The violations of scaling with mass are a consequence of the different 
possible scattering processes which lead to differing electromagnetic showering profiles. The 
tightest bounds for shorter lifetimes come from COMPTEL observations of the low energy 
part of the spectrum, while for longer lifetimes the highest observable energy at EGRET 
gives the tightest bounds. We discuss the implications of the bounds for dark matter 
candidates as well as relics that have a mass density substantially below the critical density. 
These bounds can be utilized to eliminate models that contain relics with lifetimes longer 
than 10 -4 times the age of the universe. 
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1 Introduction 



While the standard model of particle physics has passed all experimental tests to date, it 
lacks a particle candidate that could provide the dark matter of the universe as expected from 
astronomical observations. Furthermore, our present understanding of structure formation 
seems to indicate that some fraction of the dark matter should be "cold" , so as to generate the 
proper power spectrum. Such dark matter candidates are quite common in many extensions 
of the standard model. Indeed, many models predict long lived relics that may or may not 
be dark matter candidates. Long lifetimes for heavy relics, where by "long" we mean within 
several orders of magnitude of the age of the universe may arise in many models which have 
symmetries that are only broken at short distances. Thus it is interesting to investigate the 
observational signatures of such long lived relics in an effort to rule out classes of models. 

In this paper we study the signatures of particles with lifetimes comparable to the age 
of the universe. Such particles could play a role in solving the dark matter problem, but we 
will not confine our analysis to dark matter candidates. The inclusion of long lived heavy 
(M x J> 50 GeV) particles necessitates an extension of the standard model. These relics could 
be technibaryons in technicolor models or the lightest supersymmetric partner in an R-parity 
violating supersymmetric extension of the standard model. The bounds found here are model 
independent and depend on only three parameters: the mass M x , the lifetime t x and the 
radiative or hadronic branching ratio times the relic density B ■ r] x , with r] x = n x /n^. Given 
any model, it is possible to calculate the relic density using standard techniques leading to 
bounds on couplings as well as masses. 

Our exclusion bounds are derived by considering the direct observation of the gamma rays 
produced in the decay process. In general, the predicted observed spectrum will differ greatly 
from the decay spectrum due to redshifting and scattering in the early universe. The final 
spectrum can be compared to the EGRET and COMPTEL data leading to the exclusion plots 
presented here. 

Previous investigations |l|, ^, |3|, ||, |5| of the gamma ray spectra produced by long lived relics 
have concentrated on radiative decays into either photons or charged particles. We consider 
both radiative and hadronic decays^], including the effects of photon-photon scattering and 
e + e~ pair production. There are reasons to believe that hadronic decays are more compelling. 
First, the hadronic branching ratio of the relics is expected to be of order one, unless there is a 
symmetry which forbids such decays^]. Radiative decays usually arise at the one loop level and 
3 Some estimates for hadronic decays were discussed in Ref. jj]]. 

4 One would expect that the branching ratio into charged leptons may also be of order one, but then photons 
would only be generated if the lepton is energetic enough to shower. See above and the footnote on p. 423 of 
Ref. [^] . Bounds from measurements of the galactic positron flux were considered in Ref. || . 
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so one would expect smaller branching ratios. Second, hadronic decays produce more photons 
in the softer part of the spectrum due to fragmentation, which produces a large number of 
pions that decay to two photons. Thus, as opposed to the case of radiative decay, we expect 
the non-scattered spectra to have more photons at smaller energies (i.e. £L -C M x ). Since 
the diffuse photon background spectrum is well measured only up to 10 GeV, we would naively 
expect our bounds to apply to larger mass relics M x > O(10 3 ) GeV for hadronic decays but not 
for radiative decays. However, bounds derived from radiative decays of large mass relics can be 
obtained if a reprocessing mechanism to lower the photon energy is active. Such mechanisms 
begin operating once injected photon energies are above about 23 GeV where radiative decays 
can potentially compete with hadronic decays and produce large numbers of photons at low 
energy for masses larger than roughly 50 GeV. Thus, a complete calculation is necessary to 
determine the bound for a given mode of decay, as we present here. 

It should be noted that some of the bounds derived here will overlap those coming from 
structure formation arguments in the part of parameter space where the relic density is of the 
order of the critical density. Radiative decay can lead to a radiation dominated epoch after 
recombination that would drastically distort the observed power spectrum. Thus the shorter 
lifetimes considered here could lead to such distortions, but we will not consider these effects. 



2 Electromagnetic Cascades in the Early Universe 



A high energy photon injected in the early universe will, in general, scatter. Since there are 
many scattering processes possible, the nature of the scattering is strongly dependent on the 
redshift at which the photon is injected as well as its energy. Each process has a characteristic 
optical depth which determines its relevance to the evolution of the photon. The relevant 
processes were investigated in detail by Zdziarski and Svensson in Ref. Q. The processes 
include: e + e~ pair production (PP), photon-photon scattering (GG), Compton scattering and 
pair production off of matter (PPM)0. Figure |l| (derived from Q) divides the redshift energy 
plane into regions labeled by the process that dominates in a particular part of the graph. The 
photons always have injected energies, E~ > 100 MeV, which is obviously true of radiative 
decays of heavy relics and is also true of hadronic decays due to a cutoff in the spectrum at 
E~ = m n /2. Therefore, the mechanisms for rescattering photons with energies below about 
100 MeV are irrelevant to our analysis. Details of the spectra and cutoffs are described in the 
following sections. 

Given the initial energy of the photon and the redshift at which it was injected, the progress 
of the injected photon can be tracked by moving horizontally across Fig. p]. For large enough 
5 Pion production off of matter may also be of relevance for certain epochs, see Ref. H. 
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photon energies, the time scale for scattering is short compared to the expansion rate of the 
universe, so we may neglect any vertical motion until we reach the region where Compton 
scattering dominates (not shown; to the left of the left edge of the graph) or where the photon 
reaches the point where its optical depth drops below one. In the limit that the photon energy 
is much smaller than the mass of the electron (Thompson limit), the cross section for Compton 
scattering is independent of energy. When the photons reach the region of the plot dominated by 
Compton scattering in the Thompson limit, the photons will eventually be in kinetic equilibrium 
with the thermal bath leading to a finite chemical potential for the photons. The resulting 
distortion of the precisely measured microwave background leads to bounds on the injection of 
low-energy photons prior to recombination (z — 1300), as discussed in Ref. 0]. For photons with 
redshifts in the range 10 3 <> z <> 10 6 , the electromagnetic cascades result in the production 
of 3 He and D via the disintegration of 4 He. Thus, there will be even more stringent bounds for 
this range of redshifts coming from limits on primordial abundances of light elements |7|, |8L |3| . 
Here we concern ourselves with bounds coming from direct observations of gamma rays, and 
therefore we will be interested in photons whose life begins in the region dominated by pair 
production or photon-photon scattering after the epoch of recombination. Such photons will 
eventually reach the region of the plot where the optical depth drops below one and can be 
directly detected. 

Pair production leads to a cascade, since the hot electron-positron pairs produced will 
inverse Compton scatter in the Klein-Nishina regime (cr oc ln(s)/s), where they generate very 
energetic photons, which will in turn pair produce again. This process will continue until the 
photon energies drop to the point where pair production off the Wien tail of the black-body 
distribution is no longer efficient, or until we reach a regime where some other process (e.g. 
photon-photon scattering) begins to dominate. Note that during this cascade process the 
mean free path for pair production is much smaller than the local Hubble length. We may 
therefore calculate the cascade rate without regard to the effects of the expansion rate on the 
energies or occupation numbers. 

For the region of energies and redshift relevant to our analysis, photons will follow one of 
four paths. If the photon is injected in the region where the optical depth is less than one, 
then the redshifted photon will free stream to the detector. A photon injected in the region 
dominated by pair production will induce an electromagnetic cascade resulting in an "escape 



photon spectrum", that was calculated in Ref. ||, 10 1 and is presented in the next section. The 
spectrum terminates at E max , defined as the energy above which a photon will pair produce 
and degrade its energy. The line dividing the PP and r < 1 regions in Fig. [I] is E ma _ x as a 
function of redshift. A photon born in the region dominated by pair production with a redshift 
in the range 300 < z < 700 will have a slightly different fata^l. In this regime, photon-photon 



6 When 700 < z < 1300 the photons undergo pair production off of matter, but the resulting spectrum has 
not been calculated. Most of these photons will eventually reach kinetic equilibrium via Compton scattering. 
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scattering becomes relevant |J since the energy of the photons are degraded and pass through 
the region GG in Fig [I]. The width of this shaded region is given by § 

^max ^ | ^ 

where fio.i = f2&/0.1, /150 = Hq/(50 km s _1 Mpc -1 ), and T2.7 is the temperature of the mi- 
crowave background in units of 2.7 K. The existence of this region has the effect of further 
distorting the photon spectrum, as will be discussed in more detail below. Finally, if a photon 
is injected directly into the region where photon-photon scattering dominates, it will lead to 
yet another spectrum of escape photons. 




3 Scattering processes 



To quantify the photon scattering processes, we divide the region into three segments in 
redshift < z < 300, 300 < z < 700, 700 < z. This division in redshift, along with the 
forthcoming divisions in energy (E mSLX (z), -E C rit(-z)), provide an approximation to Fig. |l] that we 
use throughout the following discussion of scattering processes. As anticipated, we divide the 
region < z < 300 along the line of optical depth r = 1 for pair production, defined by || 

„ mi 36 TeV . 

^ = 30T^ TT7 (0<z< 300). (2) 
Photons injected with redshifts in the region < z < 300 with energies E~ < E max do not 
scatter, while those with energies E 7 > E max pair produce and generate a cascade. The 
resulting spectrum is given by || |9| 



0.767 E-°jE~ 1 - 5 < E 7 < 0.04£ max 

0.292 E-^E- 1 - 8 0.04£ max < E y < E max (3) 



-C/t.nt. 



-Emax < Ey 



where C is the number of photons per unit energy in the spectrum and Emax is given by (|2|). 
The spectrum is normalized according to / C(Ey)E-ydEy = E tot , where E tot is the fraction of 
energy in the injected spectrum above E max . 

In the region 300 < z < 700, photon-photon scattering dominates for photon energies 
above -E cr it and below £ max , which is given by 

Ecrit ^min 1 / , \ 

~E =- 3- ( 4 ) 

-Emax is now determined by equating the optical depths for photon-photon scattering and pair 
production. i? max is thus slightly larger || 

m 2 50 TeV 

Emax = || * ^ff^- (300 < z < 700) (5) 
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than in (||). Photons with EL < E crrt do not scatter, while those with E cr j t < E 7 < E max 
photon-photon scatter (with the background radiation), and those with E max < E 7 pair pro- 
duce. In the energy window E crit < E 7 < E max , each scattering of an energetic photon with 
a background photon results in the production of two photons which approximately share the 
energy of the injected photon. The spectrum from these scattered photons takes the form II 



2.08 



E tot E^Jl + (iW£ inj )3]V3[l + (E 7 /E crit ) 3 ] 5 /3 



(E cr it < Einj < E max ), (6) 



where E- m j is the energy of the assumed monoenergetic injected photons and Etot is the total 
integrated energy in the injected spectrum. The spectrum is normalized as in (||), where we 
assume E cr it/E max = 1/3 to obtain the overall normalization constant. A non-monoenergetic 
injected spectrum can be treated by simply splitting injected spectra into many small subregions 
between E cr i t and E max , dividing the total integrated energy in the spectrum accordingly and 
using (^). The limiting behavior of the resulting spectrum is proportional to a constant for 
E 7 <C E cr it and proportional to E~ 

spectral hump near E cr ; t . For photons injected with energy E\ 

initially scatters the photons as in (^), but photons with energy below E max can also rescatter 
by photon-photon scattering as above. The scattered spectrum can be approximated by 1|] 



for « £7 < E; n j, thus resulting in a startling 

mj > E max , pair production 



£(E 7 



E 



tot 



0.535 



10 



p2 
^crit 







H 


El f 1 




-Ecrit / 



+ 



1 + 



\ E cr it / 



5/3 



E-v 



E„ 



0.2 



0.767 E m ° x 5 E- L5 < E 7 < 0.04E, 
0-292 E-^E" 1 - 8 



0.04E n 



< E^ < E n 



(7) 



valid for E 7 < E max . 



The limiting behavior of the spectrum recovers the pair production result 
proportional to E~ L5 for E 7 < 0.04E max , while the first term in (0) dominates between 



max ^ E 7 < 0.86E n 



0.07E, 



spectrum proportional to E 



and the second term dominates for E 7 > 0.86E n 



leading to a 



In Fig. U we show the spectra for pair production in the low and high z regimes and the 
spectra for photon-photon scattering, assuming a total integrated energy of E to t in each case. 
The energy is normalized with E max = 1, E cr i t = E max /3 and E- m ^ = 2E max /3, the latter being 
an arbitrary choice (within the allowed range E max /3 < Ei n j < E max ) for illustration. 



4 Redshifting 



Given a photon spectrum after decay (from direct and/or reprocessed photons) £(E 7 ), the 
spectrum we see today is simply an integral over all redshifts convoluted with exponential decay 
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rate, 



dJ 



dE. 



70 



3_to_ 

8tt t x 



dz 



o (1 + z) 3 / 2 



exp 



A, 



1 + z)- 3 / 2 



(8) 



where 



jgjjr- is the flux of photons, E 10 is the present-day photon energy, to is the age of the 
universe, t x is the lifetime of the relic, n 7 (io) is the present-day density of photons, and n 7 
and n x are respectively the densities of the photons and relics at decoupling (we use Q, = 1 
from here on). Prom the previous section, we use zq = 700 as the upper limit in redshift since 
high energy photons injected above this redshift will pair produce off of matter and eventually 
Compton scatter thus avoiding direct detection. The photon spectrum today in (||) can be 
written suggestively as 



dJ 



dE. 



B^r] x \ 1.5 x 10 



9 z o 



To 



700 



dz 



C(E 1 



T X J (1 + 3)3/2 Br) 



r x 



,-3/2 



•cxi) ( .1 + : ) "■ ~ ^> (cm 2 s sr MeV) 

(9) 

where JC(E^) is in [GeV] -1 . It is clear from (0) that an upper limit on the present-day flux of 
photons can be translated into an upper limit on 'the ratio of relic to photon number density' 
7] x (times the relic branching ratio to photons B^) at z = 700. 

If the lifetime is short relative to the age of the universe, then most of the photons would 
be reprocessed at or before recombination, and would not reach the detector. Thus, to ensure 
that an appreciable number of photons can be observed today, we require 



r x /t Z (1 + z ) 



-3/2 



5 x 10 



-5 



(10) 



Furthermore, we have specified that the bound we obtain is the relic density at zq. This is 
roughly equivalent to the relic density at decoupling if the lifetime is longer than (|n]), so that 
the density does not change appreciably between decoupling and Zq. 



5 Diffuse Photon Background 



To establish bounds on the relic density, we use the recent bounds on the extra-galactic 
diffuse gamma ray background from the EGRET [11] and COMPTEL [12] instruments aboard 
the Compton Gamma Ray Observatory. Both instruments find that the diffuse photon flux 
obeys a power law 



dJ 



dE, 



k 



70 



R 



70 



1 MeV 



(cm 2 s sr MeV)" 



(11) 



roughly consistent with (but more sensitive than) measurements done by the SAS-2 experiment 
many years ago 13]. EGRET fit to a power law for photons in its observable range 30 < E^ < 
10 4 MeV and found k = 2.26 x 10~ 3 and a = 2.07. COMPTEL also found a reasonable fit 
with a ~ 2.3 (although they did not explicitly give fit parameters with errors). We take the 
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EGRET fit to be valid down to E J0 = 30 MeV, then we estimated a best fit power law for 
the COMPTEL data that is continuous through £ 7o = 30 MeV. We obtained a = 2.38 with 
k = 6.4 x 10~ 3 which fits the COMPTEL data quite well down to E l0 = 0.8 MeV (the lowest 
energy reported), and also fits other data ]14|] below the sensitivity of COMPTEL to roughly 
E^ =0.1 MeV. We also note here that the infamous "MeV bump" discussed in Ref. [fl has 
disappeared. 



6 Bounds from Radiative Decays 



6.1 Preamble and Previous Results 



Bounds on relics with radiative decays from diffuse background measurements have been 
considered previously in Refs. J5|. Our analysis differs in several ways. Many of the 

bounds derived here were found using the new EGRET data which allows us to look at higher 
energy gamma rays. Furthermore, our analysis of the showering profiles differs from those 
given in |[. The authors of |4], |[ determined E max by equating the Compton scattering 
cross section with that from pair production, leading to a much lower value of E max than what 
was used in this paper. This lower value of E ma _ x leads to showering at lower values of the relic 
mass and thus the bounds found in ||], |5| have a different mass dependence than found here. As 
discussed in the previous section, Compton scattering does not become important until much 
lower energies ||] and does not play a role in the determination of E may[ . We also included a 
more complete analysis of the spectral distortion due to photon-photon scattering than was 
considered in Q, though the effects on the bounds are minor. 

In what follows, we assume 2-body decays X — > 77 with a branching ratio -B 7 , giving a 
(non-scattered) input spectrum 



-C/ini hi\ ■ 



8(E^-E ini ), (12) 



^inj -^mj 

where we use E m ^ = M x /2 to represent the injected energy per decay, with a total energy in the 
input spectrum of £t t = 2-B 7 i?i n j (the factor of 2 due to two photons in the final state). The 
bounds we find for 2-body decays can be applied approximately for 3-body decays by setting 
-E; n j = (E) = M x /3. In addition, 2-body decays to single photons X — » Yj can be similarly 
constrained by scaling up the bound on r\ x by a factor of two. 
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6.2 No scattering 



In the special case that the relic decays with no scattering (so that ( |12D is the input 
spectrum), the present-day photon spectrum can be calculated exactly from (||) yielding 



dJ 



dE. 



— B 7 r/ X n 7 (i c 



To 



8tt 



3/2 



cxp 




(13) 



1/2 

for _Ei n j/(l + zq) < E l0 < E m y The photon flux rises proportional to E 10 up to roughly 
E l0 w £i n j(Tx/io) 2 ' 3 and then drops exponentially for higher i^ 7o up to E- m j. This can be seen 
in Fig. || where sample spectra are shown for E m j = 20 GeV with relic lifetimes in the range 
10 -4 < r x /t < 1. The lower bound E J0 = E- m] /(l + z Q ) ~ 28.5 MeV is clearly visible as a cutoff 
in the spectrum, as is the upper bound from the photon injection energy E 10 = E- m j = 20 GeV. 
Notice that for short lifetimes T x /to 
the final spectra for all E^ . 



10 , the exponential suppression completely dominates 



6.3 Numerical results, with scattering 



The sample spectra in Fig. ^ illustrate the effect of redshifting, but the more general case 
with scattering is what is of interest to determine relic density times branching ratio bounds. 
In Table |] we list the relevant scattering mechanisms with their redshift and injected energy 
dependence. For a given injected photon energy, we expect a dip in the spectrum due to the 
transition from scattered to unscattered photons as Ei n j is increased. The dip is located at 



Ed 



ip 



where 



and 



E clit (z=0) 



E CT it(z = 0) 



E max {z = 0) 



(£ inj < 55 GeV) 

(55 < £ inj < 122 GeV) 

(E in] > 122 GeV) 



in, 



3 • 22T 



in. 



30T 



17 TeV 



36 TeV. 



(14) 



(15) 



(16) 



The dependence on E in ^ comes from the fact that at larger energies it is possible to scatter at 
smaller redshifts. This is readily seen in the E\ n ^ = 25, 50 GeV figures of Fig. ||. In Table || 
we evaluate ( |14| ) for the injected energies that we have considered. The photons with energies 
smaller than (or to the left of) -Edip ar e those that were reprocessed by scattering, and thus 
decayed at an earlier redshift then those with energy larger than (or to the right of) £'dip> 
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Scattering mechanism(s) 






< 23 GeV 


no scattering 




23 


-> 55 GeV 


photon-photon scattering 




00 


— ► /i >je v 


photon-photon scattering 


r(l + Z> 301) 


71 


-> 167 GeV 


photon-photon scattering 


5 (301 < l + z< 50 / eV ) 


71 


-> 167 GeV 


pair production (1 + z > 


50 TeV^ 


122 


-> 167 GeV 


pair production ( 36 g TeV < 1 + z < 301) 


167 


-> 3.7 x 10 4 GeV 
> 3.7 x 10 4 GeV 


pair production (1 + z > 
pair production (all z) 


36 TeV \ 

£inj > 



Table 1: Scattering mechanisms for different injected photon energies E m j and redshifts z. 

that were unprocessed. As we further increase pair production turns on as seen in the 
Sinj = 100, 200 GeV figures of Fig. ||[ In these figures we see the power law behavior expected 
at low energies coming from those photons produced in the pair production cascade. 

Notice in particular that E- m ^ = 100 GeV displays a scattered spectra that is consistent 
with @, whereas the scattered spectra for E- m j > 200 GeV is consistent with @ (see Fig. |2|). 
The reason for this difference is that the threshold for pair production at low z(< 300) is crossed 
once Einj > 200 GeV, and consequently those scattered photons dominated the final spectra. 
It is also important to notice the critical injected energy E^ = 600 GeV corresponds to an 
-^dip = 10 GeV in Table ^. For injection energies above this value, the only detectable photons 
will be those which undergo photon-photon scattering. 

In Fig. |, we have sliced the previous photon flux vs. photon energy plots along the energy 
axis for a particular set of observed energies E l0 = 1, 10, 100, 1 x 10 3 MeV. The bound on the 
relic density can be found by using the observational limit on the diffuse background found in 
Sec. |5| for each photon energy E^ , and effectively inverting the graphs in Fig. || to give Fig. |(| 
These figures demonstrate that the best bound is not a trivial function of the measured photon 
energy, relic mass or lifetime. For example, in the E^ = 100 MeV figure one finds a better 
bound on a 25 GeV relic particle than for somewhat heavier relics. This is due to the fact 
that the unscattered photons will populate the higher energy range of the observed spectrum 
which is more strongly constrained. Whether or not such an inversion comes about depends 
upon whether the energy we are considering is larger or smaller than i?dip> defined in (|l4|). For 
shorter lifetimes and larger photon energies, one finds the unscattered part of the spectra is 
exponentially suppressed, as can be seen in the lack of a limit for EL = 1 x 10 3 MeV and small 
injected energies < 50 GeV. 
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£ inj (GeV) E dip (GeV) 



25 


0.037 


•50 


0.15 


1 nn 


U.oo 


200 


1.1 


400 


4.4 


600 


10 


800 


17 


1600 


70 


3200 


280 


6400 


1100 



Table 2: The location of the dip in redshifted spectra -Edip; where above (below) this value the 
photons originated from unscattered (scattered) injected photons (see fll4})). The entry in bold 
-^dip = 10 GeV, corresponding to an injected energy E- m ^ = 600 GeV is a critical point where 
injected energies above this value can only be detected today through scattered photons. 

Finding the maximum photon flux above background (i.e. one point on each line in each 
graph of Fig. Q), allows one to derive the bound on the relic density for a given mass and 
lifetime as shown in the upper graph of Fig. 0. Bounds for lifetimes less than T x /to ~ 10~ 5 
become poor very quickly due to the exponential suppression. Bounds for lifetimes longer than 
T x /to = 1 scale by a factor to/T x (outside (||)) relative to the bounds at r x /to = 1. We used 
the observational diffuse background fit as described in Sec. || and thus a bound for any given 
mass and lifetime utilizes one (optimal) observational energy. This is shown in the lower graph 
of Fig. [7|. For example, for r x /to > 1 one can see the trend in increasing EL, is to increase 
the photon flux (see Fig. ||). Hence the best bound for this lifetime comes from observations 
of the most energetic photons. On the other hand, for T x jt§ = 10~ 4 one finds the best bound 
for Ei n j = 25 GeV is roughly E yo 80 MeV. Higher E J0 simply pushes into the exponential 
suppression regime where no bound exists. The upper limit on M x r\ x ~ 2.5 x 10 -8 GeV 
corresponds to the critical density Q x h 2 ~ 1, which is the upper limit for any relic based upon 
the age of the universe. 

The general behavior in Fig. [/] is an increasing upper bound on M x rj x as the injected energy 
is raised up to about Ei a j = 800 GeV, and then a steady decrease thereafter for larger injected 
energies. The reason for this trend in the bounds for 2-body decays is due to the transition 
noted in Table |2] when the injected energy crosses E- m j = 600 GeV. As remarked above, the 
value -Edip = 10 GeV at this transition implies that for all injected energies above 600 GeV, 
bounds can only be derived using the scattered photons. Since the number of photons increases 
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as the injected energy increases in the scattering regime, the bounds are stronger as the energy 
is increased above 600 GeV. This effect is further enhanced by the fact that as the mass of the 
relic increases showering can occur at smaller redshifts. As the injected energy is reduced below 
600 GeV, more of the non-scattered, redshifted photons appear, and so a better bound comes 
from lowering the injected energy. This is clear since the best bound comes from the lowest 
injected energy considered, E^ = 25 GeV. The one special case is for E- m j = 800 GeV where 
the best bound for long lifetimes T x /to > 1 is not the highest present-day detection energy, 
but instead slightly less E 10 « 4 X 10 3 MeV. The reason the bound comes from lower energies is 
due the E~ 4 ' 8 suppression in the scattered photons that exists for present-day energies within 
a factor of E m3iX / E cri t = 3 lower than E^p = 17 GeV. Since the diffuse background scales as 
£>-2.07 to -2.38^ best bound will be determined by the point where the redshifted spectrums' 
slope is equal to the slope of the diffuse background . This point is roughly given by Enp/3. 

It is also interesting to note that the bounds for shorter lifetimes show scalings with the 
mass. This can be seen by the fact that for shorter lifetimes the bounds begin to lay on top of 
each other in Fig. This scaling is seen to group into two lines: one line of which consists of 
the relics whose bounds come from scattering (i.e. E^ > 600 GeV) and another line for those 
relics whose photons can be directly detected. 



7 Bounds from Hadronic Decays 



The bounds on radiative decays calculated above are strong, but it is not obvious that 
such decays ought to dominate the branching ratio of heavy relics. Here we establish bounds 
on relic particles that decay through hadronic channels. We consider 3-body decays of relics 
into all kinematically available quark pairs and one uncolored (assumed massless) spectator. 



7.1 The photon spectrum from hadronic decays 



We assume general vector-axial couplings leading to both charged and neutral current 
mediated decays. (CKM mixing is ignored since it is in general model dependent and can 
be absorbed into the couplings). We take mt = 175 GeV, and so two thresholds exist with 
increasing mass M x of the relic: M x = mt + m& (for charged current mediated decays) and 
M x = 2mt (for neutral current mediated decays). Thus, for a given mass M x , an ensemble 
of relic decays with final quark energy and momenta spanning the 3-body phase space can be 
constructed. 



The quark pairs are fragmented and decayed according to the string fragmentation scheme 15] 
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implemented in Jetset [16|. This fragmentation scheme has been well tested with collider ex- 



periment data, and the resulting photon spectrum from Jetset is discussed in Ref. [17]. In 



particular, the exact of shape and normalization of the photon spectrum depends on the par- 



ticular final state quarks [17] but generally scales with the relic mass M x , which we discuss 
below. Once the spectrum is calculated, the present-day photon flux can be determined from 

The 3-body decay allows energy to escape with the uncolored decay product. Hence, the 
integrated photon spectrum appearing after the hadronization and decay of the qq^ system 
is always less than M x . Furthermore, the hadronization process does not uniquely end with 
neutral pions, as some energy leaks into leptons. We consider only the effects of pion decays 
into photons. Therefore, the fraction of energy appearing in the final state (before redshifting) 
is generally between about M x /10 and M x /3 on average. This is to be contrasted with the 
2-body radiative decays, where the total photonic energy injected is equal to the mass of the 
relic. 



7.2 Generating the photon spectrum 



The photon spectrum is obtained directly from the hadronization of different final state 
quark pairs and is presented in Fig. ||. As we discuss below, charged current and neutral 
current mediated decays yield nearly identical results (except near the thresholds associated 
with producing one or two top quarks). As such, we will only consider neutral current decays. In 
total, 10 4 events were generated for each final state quark pair, with each event corresponding 
to a point in the 3-body decay phase space[|. We have included final state electromagnetic 
and QCD radiation showers prior to fragmentation, although the showering is performed only 
off the final state quark pair. The particular quark pair it is a special case since the top 
quark decays before it fragments. However, the photon spectrum is not particularly sensitive 
to Monte Carlo ordering of decay vs. fragmentation (Jetset fragments before decay). In 
particular, we compared the photon spectrum produced from the Pythia process e + e~ — ► it 
after fragmentation (with no initial state radiation) at large energies > 2wi( with the top 
quark decaying before and after fragmentation. The photon spectrum is slightly enhanced for 
energies E^/y/s > 0.05 when the top quark is allowed to decay before it fragments. However, 
the dominant effect on our bounds from relic decays into heavy top quarks comes from the large 
number of photons produced at low energies. In fact, Fig. || clearly shows a much larger number 
of photons from the hadronization of light quarks (dd, uu, etc.) at energies E 1 /M x > 0.01 
where one would roughly expect the slight enhancement in the top quark photon spectrum. 

7 Since the energy of the qq system has the typical 3-body distribution, a direct comparison cannot be done 
between our results and those of an annihilation signal with fixed energy as done in Ref. JT^]. However, we have 
checked that in the appropriate limit we recover the shape and normalization found there. 
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Thus in regimes where the higher energy photons define the bound, it is the light quarks' 
photon spectrum that is crucial. 

The photon spectrum originates almost entirely from decaying pions, which are created 
in the quark fragmentation process |l5|, |T^| . Thus, for a given decay M x — ► qq, the spectra 
scale with M x and can be normalized to the mass of the relic. It is only the finite mass of the 
pion that breaks the scaling behavior. The photon spectrum is also virtually independent of 
the vector-axial couplings in the 3-body phase space. The only energy dependence enters in 
the mass of the uncolored product which we take to be negligible compared with the mass of 
the relic. In fact, a large non-zero uncolored product mass could be easily accommodated by 
simply reducing the relic mass by approximately the mass of the uncolored product. 

Once the photon spectrum has been obtained from Jetset, it is fit to a sum of exponen- 
tials 0] 

C{E~) x M x = ^2 = Ae~ ax + Be^ x , (17) 
ax 

where x = E~/M x , and A, B, a, (3 are positive constants. The fit gives a reasonable character- 
ization of the photon spectrum valid for £L > m n /2. The dependence of the fit on the Monte 
Carlo statistics is small. Increasing the statistics by a factor of 5 shifts the final redshifted 
spectra by at most 0(15%). 

7.3 Numerical results: Photon flux versus photon energy 

We have scanned the relic mass range over more than two orders of magnitude from M x = 
50 GeV (a likely lower bound from LEP) through 12.8 TeV in steps of a factor of 2. We have 
simultaneously scanned the lifetime range throughout the region that gives bounds for relevant 
densities. Fig. || displays a selection of the above sampling, with relic masses M x = 50, 100, 
800, 6400 GeV and lifetimes including r x /t = 10" 4 , 10" 3 , 10~ 2 , 10 _1 , 1. The figures show 
the photon flux as a function of the present-day photon energy. We have also calculated the 
photon flux as a function of the photon energy for the same relic mass-lifetime parameter space 
using charged current interactions. The photon fluxes are virtually identical throughout most 
of the mass range, with the neutral current interaction usually giving a slightly larger value 
(due to decays into top quark pairs) than decays via charged current interactions. However, in 
the mass window mt + mj, < M x < 2 nit the photon flux from a charged current interaction 
is about a factor of two greater than the equivalent spectra from a neutral current interaction, 
since the bt decay mode is open. 

It is clear from Fig. || that most of the photons are well below M x /10, so that the effect 
of scattering for hadronic decays is not as prominent as in radiative decays for the same relic 
mass. The effects of photon-photon scattering and pair production are handled similarly to 
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the radiative decays, however the input spectra is no longer a delta function in the photon 
energy. Specifically, the spectra for photon-photon scattering @, pair production at low z @ 
and pair production at high z (0) are used in the same way as in radiative decays. We need 
only determine how much energy is injected into each regime. Since the energy injected into 
this regime is dependent on -E max and E cr n, which are in turn dependent on the redshift z, the 
procedure must be done numerically. 

As noted above the non-scattered spectra are cutoff at E^ = m n /2, whereas the scattered 
spectra have no such cutoff. Thus, in regimes where the non-scattered (injected) spectra dom- 
inate the final redshifted spectra, we expect a cutoff at E l0 = ^/(l + zq) ps lO^ 1 MeV. Such 
a cutoff is observed in all of the spectra in Fig. [9|. Fig. || also shows explicit scaling with relic 
mass in the redshifted spectra. This is a consequence of the string fragmentation process, which 
is Lorentz invariant (given that enough energy is present to create a string and subsequently, 
jets). Scaling does not hold, however, for the scattered spectra since there are absolute cutoffs 
(i? max , E cr it) involved. Note also that despite the fact that there is no scattering, the location 
of the peak photon flux increases above 10 -1 MeV as the lifetime is increased. This is to be 
expected since as the lifetime T x /to — ► 1, the peak flux should approach m n /2 since photons 
that have not been appreciably redshifted (z ~ 1) are not exponentially suppressed. 

Different final state quark pairs give rise to different final spectra. However, the principle 
differences between decays into particular final state quark pairs is not difficult to understand. 
As can be anticipated from the injected (non-redshifted) spectra in Fig. ||, decays into top 
quark pairs yield the best bound when the lowest energy photons from the injected spectra are 
sampled. The decay into top quark pairs gives the largest flux of photons for relics with a large 
mass M x > 2m t at very low present-day energies E 10 <C M x /(1 + zq). Similarly, it is the decay 
into light quark pairs dd and uu that give the largest flux of photons for higher energy photons 
E~f > M x /(1 + zq). Without a theoretical motivation for decays into one or another flavor 
or family, we choose to divide the branching ratio equally among the kinematically available 
quark pairs. Thus we assign an equal branching ratio for all the pairs (1/5 or 1/6 depending 
on whether the tt threshold has been crossed). 

The effect of scattering on the spectra, as remarked near the beginning of this section, is 
not as important for hadronic decays as it is for radiative decays. In fact, scattering is virtually 
absent for M x = 50 GeV, as illustrated in Fig. [9| where there is no photon flux (above the lower 
limit in the graph) for E~ < 1CP 1 MeV. This is not surprising since the quark pair will always 
have an invariant mass less than M x /2, which is only barely above the threshold for photon- 
photon scattering ~ 23 GeV. For M x = 100 GeV, only photon-photon scattering is possible, 
and one can see the characteristic limiting behavior of a flat spectra for E~ Q < 10 _1 MeV. 
For M x > 200 GeV, pair production dominates the scattered piece of the redshifted spectra 
E l0 < 10 _1 MeV, albeit with a total integrated energy that is much less than the unscattered 
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piece (E yo > W~ 1 MeV). It is really only when M x > 1 TeV that the scattered piece of the 
spectra has a photon flux comparable to the unscattered piece. This implies that the bulk of 
the injected photons are below about M x /10, which is roughly the scale where scattering turns 
on. 



7.4 Relic density bounds 



In Fig. [H] we have sliced the previous photon flux vs. photon energy plots along the energy 
axis in analogy to Fig. ||, for the particular present-day energies E l0 = 1, 10, 100, 1 x 10 3 MeV. 
Just as in the 2-body case, the bound on the relic density can be found by using the observational 
limit on the diffuse background found in Sec. |5] for each photon energy E^ , as shown in Fig. 11. 
We observe, as in the 2-body radiative decays, that the mass dependence shows nontrivial 
behavior characteristic of the transition between relics whose bounds come the scattered and 
unscattered spectra respectively, for M x > 1 TeV. In addition, one finds that for shorter 
lifetimes and larger photon energies the spectra are exponentially suppressed, as can be seen 
in the lack of a bound for £L = 10 3 MeV and small masses M x < 100 GeV. The physical 
interpretation is that most of the decays occurred much earlier than our present epoch, so the 
photon flux is significantly more redshifted than for relics with a longer lifetime. Thus, we see 
a much smaller number of present-day photons at high energy. 

By finding the maximum photon flux above background (i.e. one point on each line in 
each graph of Fig. ||), one can derive the bound on the relic density for a given mass lifetime 
as is done in Fig. ^ (the same procedure as in Fig. |^). We used the observational diffuse 
background fit as described in Sec. ||, and thus a bound for any given mass and lifetime utilizes 



one (optimal) observational energy, as is shown in the lower graph of Fig. 12. For example, for 
T x /to = 1 one can see the trend in increasing E l0 is to increase the photon flux (see Fig. |9|). 
Hence the best bound for this lifetime comes from observations of the most energetic photons. 
On the other hand, for r x /to = 10 -4 , one finds that the best bound for M x = 50 GeV is 
roughly E 10 rj 7 MeV. 



8 Implications of the Bounds 



There are two central results we can extract from Figs. and 12 for both radiative and 
hadronic decays. First, a large range of lifetimes can be excluded for a relic with roughly 
the critical density. Second, relics with densities considerably below the critical density are 
excluded, which places a strong constraint on models with a long lived massive particle. 
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8.1 Relics with the critical density 



For relics with roughly the critical density that decay dominantly through a radiative 
channel, the bounds from the diffuse background exclude lifetimes in the range 

10 12 < t x < 3 x 10 22 s (18) 

(using to = 10 10 yr). This bound applies to a relic with any mass M x > 1 MeV. The upper 
bound on the excluded lifetime of 3 x 10 22 s applies to the worst-case scenario with M x /2 = 
-E'inj ~ 600 GeV, where the upper bound increases to roughly 10 27 s for M x /2 = £ inj = 25 GeV 
and to roughly 10 25 s for M x /2 = E[ n j = 6400 GeV. The upper bounds increase for larger mass 
relics M x J> 10 TeV. The upper bound on the excluded range is near 10 28 s for small masses 
0(1 GeV). 

For relics with roughly the critical density that decay dominantly through hadronic chan- 
nels, the bounds from the diffuse background exclude lifetimes in the range 

10 12 < r x < 10 26 s (19) 

for masses M x = 50 — > 10000 GeV. Unlike 2-body decays, the upper bound on the lifetime 
decreases as the mass of the relic is increased. This is because for masses beyond 10 TeV, an 
increasing fraction of the photons are scattered into lower energies where the diffuse photon 
bound is weaker. For masses smaller than 50 GeV the upper bound on the lifetime is roughly 
similar, but is absent once the mass is below the threshold for pion production. 



8.2 General bounds 



We have shown that a large range in relic lifetime can be excluded assuming the relic 
has roughly the critical density. However, Figs. ^ and 12 clearly show that the upper limit 
on the relic density is much smaller that the critical density by several orders of magnitude. 
In particular, models with a long lived particle that does not have the critical density will be 
excluded if the relic density exceeds our bounds. The translation of our bounds is relatively 
straightforward, if most of the relics have not decayed prior to the earliest redshift which we 
considered, zq = 700 (that is, if the lifetime is longer than that in (|To|)). Specifically, 



2 Mxju 

2.5 x 10- 8 GeV V ; 



for r x /t > 10" 
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9 Conclusions 



Utilizing the latest observations of the diffuse photon background found from EGRET 
and COMPTEL leads to the bounds summarized in Figs. [7| and |l^. Since the diffuse photon 
background is now well measured (up to 10 GeV), these bounds establish a fixed upper limit on 
the relic density of long lived relics. We find that 3-body hadronic decays typically give better 
bounds than 2-body radiative decays for the same relic mass despite a larger total energy 
deposited into the final spectra for the latter decay mode. The stronger bounds on hadronic 
decays are a direct result of the low energy photons emitted from the fragmentation process that 
produces pions which then decay into two photons. However, strong limits on radiative decays 
have also been obtained for heavy mass relics since high energy photons are typically scattered 
by either photon-photon scattering or pair production. In particular, relics with the critical 
density and the masses considered here that decay dominantly through radiative (hadronic) 
channels are excluded for lifetimes in the range 10 12 < r x < 10 22 (10 26 ) s. The upper bound 
on the excluded lifetime assumes the worst-case, which is not necessarily the smallest or largest 
mass. For particular masses or lifetimes, considerably more stringent bounds on the relic density 
can be read off from Figs. 0, [l2|. 

The existence of strong bounds for both radiative and hadronic decays from the diffuse 
photon background provides a useful tool for those who consider long lived relics in particular 
particle physics or cosmological models. If we assume that the relic has roughly the critical 
density, then we have seen that the lifetime must be far greater than the age of the universe. 

The bounds derived here allow for a more general analysis in that we do not make any 
assumptions about the number density of the relic. Thus, whether or not the relic is the dark 
matter, one is still able to put strong constraints on the model. Furthermore, it may be the 
case that the model under consideration has more than one (meta-)stable relic, one of which 
is a dark matter candidate. Such scenarios may arise in cases where there exist symmetries 
which are only broken by higher dimensional operators. For instance, in many supersymmetric 
models it is assumed that R parity is classically conserved. In such cases one may expect that 
this symmetry will be broken by gravitational effects leading to very long lived particles which 



may or may not be dark matter candidates. Bounds in such scenarios were discussed in Ref. [18] 
using limits from the positron flux in our galaxy for the case of critical density^. These bounds 
are useful for limiting the values of the couplings involved in the decay, which in general are 
uncorrelated to all couplings that can be measured in accelerator experiments. However, if we 
do not make any assumptions regarding the energy density of the relic then we may constrain 
couplings that are accessible in collider experiments via a calculation of the relic density. The 



8 If the R symmetry is continuous, then there is a host of other constraints arising from Majoron produc- 
tion liil. 
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relic density may be determined by calculating the annihilation cross section which will be a 
function of the accessible parameters of the theory. The drawback to this scenario is that one 
must make some assumption about what is considered to be unnaturally small for the symmetry 
breaking couplings. 
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E y [MeV] 



Figure 1: Dominant scattering mechanisms for high energy photons injected in the post- 
recombination era. The region below the solid line has an optical depth r < 1, thus no scattering 
occurs. The other regions are dominated by e + e~ pair production (PP), photon-photon scat- 
tering (GG), pair production in matter (PPM1) and pair production in ionized matter (PPM2). 




Figure 2: Scattered spectra from pair production at low z, photon-photon scattering, and pair 
production at high z. Each spectra has a total integrated energy of E tot , with E max set to 
unity. See the text for details. 
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Figure 3: Photon flux as a function of photon energy (z = 0) from a non-scattered 2-body 
decay (see (12)) with M x /2 = E- m \ = 20 GeV. The spectra scale linearly with the radiative 
branching ratio £> 7 and the relic density r] x , thus an observational limit on the photon flux can 
be translated into a limit on B^rj x . A sample of relic lifetimes r x /to were chosen and plotted 
as separate contours on the graph. 
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Figure 4: Redshifted photon spectra for a relic with mass M x /2 = E{ n j = 25, 50, 100, 200 GeV, 
a 2-body decay X — > 77 with branching ratio 2? 7 , and a relic density rj x - A sample of relic 
lifetimes T x /to were chosen and plotted as separate contours on the graph. 
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Figure 5: Slices of Fig. || with present-day detection energies E JQ = 1, 10, 100, 1 x 10 3 MeV. To 
focus on the behavior of the photon flux for different masses and present-day photon detection 
energies, we restricted the lifetime r x /to to be in the range 10~ 4 — ► 1. 
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(a) E = 1 MeV (b) E = 1 MeV 
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Figure 6: Upper bounds on the relic density (times the radiative branching ratio) for particular 
present-day detection energies E 10 = 1, 10, 100, 1 x 10 3 MeV. For a given relic mass, the region 
with relic density larger (or above) the mass contour is excluded. Notice that better bounds do 
not necessarily come from higher or lower detection energies. 
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Figure 7: The upper graph shows the final relic density bound for 2-body radiative decays 
with lifetimes in the indicated range. The bound scales linearly with the radiative branching 
fraction of the relic £? 7 , although a branching ratio different from one does not strongly affect 
our bounds. The upper limit on the relic density of ~ 2 x 10~ 8 GeV is roughly the critical 
density corresponding to £l x h? ~ 1. The lower graph shows the optimal photon detection 
energy to obtain the best bound for a given lifetime. This graph is divided at E J0 = 30 MeV 
with a dotted line to show which instrument provides the diffuse photon background bound for 
a given lifetime. 
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Figure 8: The photon flux from selected final state quark pairs from a 3-body decay of a relic. 
Notice the photon spectrum for it has an excess of very low energy E y /M x <C 0.01 photons 
and a deficit of higher energy photons E^/M x > 0.01, compared with the photon spectra from 
lighter quark pairs. Note also that the spectra are effectively cutoff at E^ = m w /2, although 
the mass used to generate these spectra M x {= 800 GeV) 3> m n , and so the cutoff is not visible. 
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(a) M x = 50 GeV 



(a) M x = 100 GeV 




Figure 9: Redshifted photon spectra for a relic with mass M x = 50, 100, 800, 6400 GeV, a 
hadronic branching ratio B^, and a relic density r\ x . All kinematically allowed 3-body hadronic 
decays X — > qq (+ uncolored product) are assumed to occur with equal branching ratio. A 
sample of relic lifetimes T x /to were chosen and plotted as separate contours on the graph. 
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(a) E = 1 MeV (b) E = 1 MeV 
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Figure 10: Slices of Fig. || with present-day detection energies E l0 = 1, 10, 100, lx 10 3 MeV. To 
focus on the behavior of the photon flux for different masses and present-day photon detection 
energies, we restricted the lifetime T x /to to be in the range 10~ 4 — > 1. 
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Figure 11: Upper bounds on the relic density (times the hadronic branching ratio) for particular 
present-day detection energies E 10 = 1, 10, 100, 1 x 10 3 MeV. For a given relic mass, the region 
with relic density larger (or above) the mass contour is excluded. Notice that better bounds do 
not necessarily come from higher or lower detection energies. 
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Figure 12: The upper graph shows the final relic density bound for 3-body hadronic decays 
with lifetimes in the indicated range. The bound scales linearly with the hadronic branching 
fraction of the relic B^, although a branching ratio different from one does not strongly affect 
our bounds. The upper limit on the relic density of ~ 2 x 10~ 8 GeV is roughly the critical 
density corresponding to £l x h? ~ 1. The lower graph shows the optimal photon detection 
energy to obtain the best bound for a given lifetime. This graph is divided at E J0 = 30 MeV 
with a dotted line to show which instrument provides the diffuse photon background bound for 
a given lifetime. 
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